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ABSTRACT 


An  analysis  was  undertaken  to  better  understand  the  phenomena 
occurring  in  spray  coolers  and  to  develop  a  mathematical  model  for 
'comparison  with  experimental  data  from  an  operating  unit.  The 
physical  characteristics  of  an  operational  exhaust  gas  spray  cooler 
and  the  instrumentation  systems  are  described.  A  mathematical  model 
of  a  spray  cooler  was  developed  by  assuming  kinetic  and  thermodynamic 
equilibrium  and  one -dimensional  flow.  A  mathematical  model  of  a 
hypothetical,  optimum  cooler  is  included  in  order  to  have  a  basis  for 
defining  cooler  efficiency.  The  equations  were  programmed  for 
numerical  solution  on  a  digital  computer,  and  several  trial  case  runs 
are  presented.  Experimental  measurements  are  compared  with  the 
efficiencies  predicted  by  the  mathematical  models. 
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SECTION  I 
INTRODUCTION 


Testing  of  turbojet  engines  and  rocket  motors  at  simulated  altitude 
in  ground  test  facilities  requires  cooling  of  the  high  temperature 
exhaust  gas  to  a  relatively  low  temperature  before  the  gas  enters  the 
exhaust  gas  pumping  system.  Cooling  of  the  gases  by  water  spray  with 
direct  heat  and  mass  exchange  between  the  water  and  the  exhaust  gas 
has  been  utilized  in  many  test  facilities.  This  method  of  cooling  is 
often  called  spray  cooling. 

Many  of  the  spray  coolers  used  in  the  Engine  Test  Facility  (ETF) 
at  the  Arnold  Engineering  Development  Center  (AEDC)  receive  exhaust 
gas  from  a  rocket  or  turbojet  engine.  The  cooling  process  reduces  the 
temperature  from  approximately  4000°R  (maximum  temperature  of  a 
turbojet  engine  exhaust  gas)  to  approximately  550°R.  By  means  of  an 
atomizing  water  spray,  the  exhaust  gas  is  cooled  and  dehumidified. 

The  cooling  produces  a  temperature  compatible  with  the  ducting,  con¬ 
trol  valves,  and  pumping  system  material  limits,  and  the  dehumidifica¬ 
tion  maximizes  the  exhaust  gas  handling  capacity  of  the  exhaust  pumping 
machinery.  Water  conservation  is  an  important  consideration  in  opera¬ 
tion  because  of  the  large  quantities  of  spray  water  required. 

There  is  limited  knowledge^of  the  detailed  performance  of  exhaust 
gas  coolers  in  general  and,  in  particular,  of  coolers  operating  at  low 
pressures  and  with  high  velocity  gas  streams.  Therefore,  a  program 
was  initiated  with  two  objectives:  ( 1)  to  develop  a  mathematical  model 
to  assist  in  the  formulation  of  design  criteria  for  future  coolers,  and 
(2)  to  define  and  understand  the  general  operating  characteristics  of  the 
coolers  in  use  at  the  ETF.  When  data  collected  during  exhaust  gas 
cooler  operation  in  support  of  normal  engine  test  programs  are  com¬ 
pared  with  data  from  the  mathematical  model,  the  limits  of  operation 
for  the  coolers  currently  in  use  can  be  determined. 

The  overall  work  is  being  carried  out  as  a  joint  effort  between  the 
Arnold  Engineering  Development  Center  (AEDC)  and  the  German  Agency 
for  Aeronautics  and  Space  Research  (DFVLR).  The  work  at  AEDC  was 
conducted  in  the  Engine  Test  Facility.  The  cooler  performance  data 
were  obtained  from  one  of  the  propulsion  engine  test  cells  during  actual 
test  runs.  Instrumentation  for  the  cooler  performance  data  was  in¬ 
stalled  in  the  exhaust  gas  cooler  and  associated  ducting.  The  work  at 
DFVLR  (Refs.  1  and  2)  was  performed  at  the  Institute  for  Chemical 
Rocket  Testing  at  Lampoldshausen.  The  DFVLR  work  is  developed  for 
low  velocity,  low  mass  flow  rates  of  exhaust  gas,  whereas  the  work  in 
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ETF  is  primarily  concerned  with  high  velocity,  high  mass  flow  rates  of 
exhaust  gas. 

This  report  contains  the  derivation  of  a  mathematical  model  based 
on  thermodynamic  and  kinetic  equilibrium  at  the  cooler  exit.  A  cooler 
efficiency  based  on  a  hypothetical  optimum  cooler  was  defined.  The 
physical  characteristics  of  an  operational  exhaust  gas  spray  cooler  and 
the  instrumentation  systems  are  described.  Experimental  measure¬ 
ments  obtained  from  actual  cooler  operation  are  compared  with  the 
mathematical  model  predictions.  The  mathematical  model  was  used 
over  a  full  range  of  inlet  conditions  to  calculate  cooler  exit  conditions, 
and  the  results  are  presented.  A  bibliography  of  related  reports  is 
also  presented. 


SECTION  II 

EXHAUST  GAS  COOLER  SYSTEM 


2.1  CONFIGURATION 

The  configuration  for  a  typical  horizontal  exhaust  gas  spray  cooler 
consists  of  a  diverging  conical  inlet  section  followed  by  a  constant-area 
duct  to  the  end  of  the  spray  cooler  (Fig.  1,  Appendix  I).  The  cooling 
water  is  introduced  through  groups  of  nozzles  arranged  in  a  series  of 
banks  (Fig.  2a),  in  which  the  first  several  banks  of  sprays  consist  of 
nozzles  projecting  a  fan-type  spray  directed  downstream  along  the  wall 
to  protect  the  ducting  (Fig.  2b).  The  other  banks  consist  of  spray  heads 
arranged  in  a  wagon-wheel  configuration  with  several  spokes,  each 
"spoke"  containing  several  spray  heads.  Each  spray  head  contains 
several  fixed-geometry,  conical  spray  nozzles  (Fig.  2c)  directed 
generally  downstream.  The  water  is  supplied  by  a  large  header,  and 
each  spray  bank  is  supplied  from  the  header  by  a  line  containing  a  valve 
for  flow  control. 


2.2  INSTRUMENTATION 

Instrumentation  was  provided  to  measure  flow  rates,  pressures, 
and  temperatures  of  the  exhaust  gas  stream  entering  and  leaving  the 
cooler;  the  cooling  water  into  the  spray  cooler;  and  the  drain  water  out. 
of  the  cooler  at  the  cooler  drain  line.  The  location  of  this  instrumenta¬ 
tion  is  shown  in  Fig.  3.  Measurements  taken  by  this  instrumentation 
provided  information  which  was  used  to  provide  experimental  correla¬ 
tion  with  the  analytical  results  from  the  mathematical  model. 
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The  millivolt  outputs  from  the  thermocouples  and  strain-gage-type 
pressure  transducers  were  recorded  on  either  magnetic  tape  or  by  a, 
photographically  recording  galvanometer-type  oscillograph.  The  mag¬ 
netic  tape  data  were  reduced  on  a  digital  computer,  and  the  oscillograph 
data  were  reduced  manually  using  electrical  calibrations  taken  prior  to 
testing. 

The  cooling  water  flow  rates  were  calculated  using  the  spray  noz¬ 
zle  manufacturer's  value  of  discharge  coefficient  and  a  measured  pres¬ 
sure  drop  between  the  supply  header  around  the  cooler  and  an  average 
cooler  pressure.  This  method  assumes  that  the  discharge  coefficient 
and  the  area  of  the  nozzles  have  not  changed  appreciably  from  the 
original  values.  The  measured  pressures  were  recorded  by  an  oscillo¬ 
graph. 

The  cooling  water  leaving  the  exhaust  gas  cooler  at  the  exit  through 
the  large  (14-in.  -diam)  drain  line  was  measured  using  a  small  flow¬ 
meter  (3 /4-in.  turbine-type)  mounted  as  a  probe  in  the  drain  line.  The 
data  were  recorded  by  (1)  an  oscillograph  as  an  analog  signal  and 
(2)  a  frequency  signal.  The  performance  of  this  device  wa§  not  satis¬ 
factory  because  of  the  very  low  flow  rates  during  cooler  operation. 


SECTION  III 

DEVELOPMENT  OF  THE  ANALYTICAL  MODEL  OF  THE  SPRAY  COOLER 


The  model  to  predict  cooler  exit  conditions  assumes  that  the  exhaust 
gas,  water  vapor,  and  liquid  water  are  in  thermodynamic  equilibrium 
and  have  the  same  velocity  (kinetic  equilibrium)  at  the  exit.  The  deriva¬ 
tion  of  the  equations  necessary  to  compute  the  exit  conditions  and  effi¬ 
ciency  of  the  exhaust  gas  cooler  is  outlined  below. 


3.1  DERIVATION  OF  EQUATIONS  GOVERNING  COOLER  BEHAVIOR 

From  given  cooler  inlet  conditions  (pressure,  temperature,  gas 
velocity  and  composition  and  cooling  water  temperature  and  velocity), 
the  amount  of  cooling  water  necessary  to  reach  a  preselected  exit  tem¬ 
perature  is  computed  in  two  distinct  processes:  (1)  saturation  (evapora¬ 
tion)  and  (2)  dehumidification  (recondensation)  with  the  gas  remaining 
saturated.  The  exit  pressure,  velocity,  and  volume  flow  rate  of  non¬ 
condensable  gas  is  then  calculated.  The  division  of  the  cooler  into 
distinct  processes  is  an  idealization  of  the  real  process  and  is  done  for 
analytical  reasons  only.  A  schematic  of  the  cooler  is  shown  in  Fig.  4. 
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The  analysis  of  the  cooler  is  based  on  the  following  assumptions: 

1.  Thermodynamic  and  kinetic  equilibrium  exist  between 
exhaust  gas,  water  vapor,  and  liquid  water  at  the 
cooler  exit. 

2.  Exhaust  gas  and  water  vapor  behave  as  perfect  gases. 

3.  All  cooling  water  is  injected  in  the  axial  direction  and 
remains  in  the  stream  as  vapor  or  liquid. 

4.  The  flow  is  one  dimensional,  and  the  cooler  has  a  con¬ 
stant  cross-sectional  area. 

5.  Drag  forces  on  duct  walls  and  spray  water  piping  are 
negligible. 

6.  All  cooling  water  is  injected  at  the  same  temperature, 

Twj  =  TW2‘ 

The  necessary  equations  are  those  of  conservation  of  mass,  species, 
momentum,  and  energy  and  the  equation  of  state.  Figure  5  is  a  schematic 
drawing  showing  the  process. 


i^ncj 

+ 

mw  j  +  mv  j. 

=  n>nc2 

"•nc  j 

=  "'nc 

~  Pnc 

mw  J 

* 

rnVj  =  mv  ^ 

For  the  saturation  (evaporation)  process: 
Conservation  of  mass: 

also 
and 

Conservation  of  momentum: 

PjA  +  (mVj  +  nine  j)  ncj  +  "’wj  ^ w  i 

=  P2A  +  (niv2  +  n,nc2^  ^  ncjj 

Conservation  of  energy: 

"lvi(hvi  +  “if-)  +  (hnCl  +  +  "*Wl  (hwi  + 

.  (  W\  .  (  . 

=  mVj  ^>V2  +  — - - J  +  mnc2  ^hnc2  +  2  I 


) 


where  assumption  (1)  has  been  used. 


(1) 

(2) 

(3) 

(4) 


(5) 
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Equation  of  state: 


(6) 


For  the  saturation  process,  energy  Eq.  (5)  can  be  written  as: 


where  Eq.  (3)  has  been  used  to  eliminate  mw^ 


For  the  dehumidification  process: 
Conservation  of  mass: 


also 


mnc  2  **■  mw 2  **■  mV2  —  mncg  +  +  ^V3 

^nc2  =  ”1nC3  =  ™nc  =  Pnc  AV nc 


3w2  +  mv  2  —  mwg  +  mv3 


Conservation  of  momentum: 

P2A  +  (mV2  +  mnC2)  VnC2  +  mW2  V„2  =  P3A  +  (mV3  +  mnC3  +  mW3)  Vnc-3 

/ 

Conservation  of  energy: 


Equation  (6)  can  be  written  in  the  form: 

Pnc  2  PnC2  ^  3 

Pn  =  3  P  n  c  3  T2 

Combining  Eqs.  (2)  and  (8),  squaring,  and  solving  for  VnCg2  give 


(8) 


(9) 
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By  noting  that  m  =  pVA  and  using  assumptions  (1)  and  (4),  the  following 
can  be  written: 


v3 


v3  “‘v 


P  nc 


Pncs  mnc 


(10) 


By  solving  Eq.  (3)  for  mWg  and  substituting  Eqs.  (9)  and  (10)  along  with 
mW3  into  Eq.  (5),  the  ratio  m^/n^nc  can  be  expressed  as: 


mw2  ^nC2  ^nc3^“ 

II!)  - 

^nc2  rav2  ^ V3mv 

2  mnc  (h«3  hv2)  h*3) 

mnc 

^  w  3  “  h  w  2^ 

/Pnc2T3\2  /'Vw2>l? 

L\Pnc3T2/  \VnC2/J 

V  2 
v  nc2 

2 

(ID 


This  is  the  form  of  the  energy  equation  to  be  used  for  the  dehumidifica¬ 
tion  process. 


To  develop  a  usable  form  of  the  momentum  equation,  using 
Eqs.  (2),  (3),  and  (4)' and  dividing  by  mnc  produces  the  expression: 


P3  =  P2  +  2q 


[(-£-)( 


(12) 


where 


Pnc2  ^  nc2 


Substitution  of  Eq.  (8)  into  this  equation  gives 


P3  =  P2  +  2q  17,  +  (\  -  -^lL3V  n_ 

A  mnc  f  \  ^1103  ^2/  ranc  \^n 


*2 
c  2 


P  nc2  "**3  V 

Pnc3  T2/_ 


(13) 


The  partial  pressure  of  the  noncondensable  gases  at  the  cooler  exit 
is  eliminated  by  applying  Dalton's  law  of  partial  pressures: 


P3  =  P 


V3  T-  Pnc3 


(14) 


At  this  point  the  momentum  equation  is  rearranged  and  placed  into 
quadratic  form: 


aP3^  »  bP3  c  =  0 


where 


=  1 


(15) 


=  -  F 


P2  +  Pv3  +  2q [ 1  + 


(' 


w  2 


mw2  V 


'-)] 


(16) 
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and 


c  =  2q 


w  2 

1  +  -  + 


m  vt  2  ^  w  2 


ln  n  c  ^  n  c  2 


-  2q 


p  i' 

*  nc «  1 


2^3  /m 

T2  \  m 


The  roots  of  the  quadratic  have  the  form 


P3  = 


-b  +  s/h*  —  4ac 
2a 


P2  Pv3  { 16) 


(17) 


where  only  the  {+)  sign  has  physical  meaning.  In  summary,  the  equa¬ 
tions  to  be  solved  are  Eqs.  {7)  or  (11)  and  (17)  and  the  following  condi¬ 
tions  are  known: 

l;  The  enthalpy  of  all  constituents  as  a  function  of  tem¬ 
perature  only. 

2.  The  partial  pressure  of  water  as  a  function  of 
temperature. 

3.  The  temperature,  pressure,  and  velocity  of  the 
exhaust  gas  entering  the  cooler,  and  the  tempera¬ 
ture  and  velocity  of  the  cooling  water  being  sprayed 
into  the  gas  stream. 

The  procedure  used  by  the  computer  to  solve  the  equations  is  given  in 
Appendix  III. 


3.2  "OPTIMUM"  COOLER 

An  "ideal"  cooler  may  be  defined  by  the  same  assumptions  as  used 
for  the  equilibrium  model  except  that: 

1.  The  static  pressure  will  be  constant  in  the  cooler,  and 

2.  All  liquid  water  { recondensed  and  sprayed)  will  be  re¬ 
moved  from  the  stream  when  it  has  achieved  thermo¬ 
dynamic  equilibrium. 

A  theoretical  cooler  efficiency  may  then  be  defined  as  the  ratio  of  the 
cooling  water  required  for  this  ideal  cooler  to  that  determined  by  the 
equilibrium  calculations  at'conditions  of  equal  volume  flow  rates  at  the 
cooler  exit. 

mw  ideal 

*?theo  =  T  (  18) 

mw  equilibrium 

The  requirement  of  equal  volumetric  flow  rates  was  made  because  this 
condition  is  believed  to  agree  with  conditions  in  ETF. 
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The  equilibrium  analysis  has  been  presented  in  the  previous  section. 
The  idealized  equations  are  developed  in  this  section  and,  because  of 
their  ultimate  use  in  determining  cooler  theoretical  efficiency,  just 
described,  are  called  "efficiency  equations." 

The  idealization  assumes  that: 

1.  The  cooler  operates  with  a  constant  static  pressure 
(Pi  =  P2  =  P>*  and 

2.  All  previous  assumptions  are  valid  except  that,  after 
saturation  is  initially  achieved,  all  liquid  water  is 
removed  from  the  stream  after  coming  to  thermo¬ 
dynamic  equilibrium. 


Writing  mass  and  heat  balance  equations  for  the  cooler  gives: 
Mass  balance: 

^nc  j  mVl  mw ,  —  +  111  x 

but 


ri  f  ITIwj  _  mnc2  "*■  111  v  2 


therefore. 


mnc  j  “  mnc 2 


'l  ^  ,UWJ 


Energy  balance: 

mnc  ^nc  1  +  mv1  +  niWlh,Vj  =  mnc  hnc2  +  I"v2  ^v2 

Equations  (21)  and  (22)  can  be  combined  to  yield 

hnC2  —  hnCj  =  ftiyg/ninc  (^vj  “  ^v2)  mwj-,/mQc  (I'wj  “  ^vj) 

Substituting  Dalton's  Law  [Eq.  (14)]  into  Eq.  (10)  gives 

m  y  p  v  2  P  v  2  ^  v 


(19) 

(20) 

(21) 

(22) 

(23) 

(24) 


mnc  Pnc  2  mnc 

Using  the  continuity  equation  (21)  and  solving  for  the  ratio  P^/P  give 


v2 


(m»l  m  v l \ 

~  ■*"  ",  ) 
mnc  mno  / 

(®wj  m v  j  \ 

:  "*■  T  J  mnc 

mnc  mnc / 


(25) 
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Substituting  Eq.  (24)  into  Eq.  (23)  gives 


h"2 '  l,"e‘  ■  ip.pji 


^  Vn  ^  V  mw. 

(hvl  -  hV2)  +  (hWl  -  hVl) 


V2'  111  nc 


m 


Equation  (25)  can  be  solved  for  — 


w  i 


mnc 


{P v / m v  mvA  m v 

P  ) \mnc  “nc/  mnc 


v  2 


-1 


(26) 


(27) 


This  ratio  is  used  in  Eq.  (26),  and  PV2  is  solved  for: 


pv2  =  P 


hnc2  hncj  +  rfi 


j  “  h  vj) 


m V  /  mv  01  V  j\ 

“nc  2  ~~  ^ncj  ~  ftv|  ™  "*"  (  I  ) 

mnc  \ranc  mnc  / 


(28) 


The  dehumidification  section  for  the  ideal  cooler  is  represented  by 
a  series  of  differential  elements  in  which  water  is  allowed  to  reach 
thermodynamic  equilibrium  with  the  gas  phase  and  the  recondensed 
water  is  removed  before  advancing  to  the  next  differential  element. 
Each  differential  element  corresponds  to  a  dT.  A  schematic  diagram 
is  shown  in  Fig.  6  for  a  single  typical  element. 

For  a  differential  element,  the  equations  for  the  conservation  of 
mass  and  energy  can  be  written: 

Mass: 

^*2  dmv  +  dmWjj 

Energy: 

mnc  dhnc  —  j  dmW2  —  niv  dhv  —  hvdmv  —  hWjj  dmWQ 

Substituting  for  dmWD  from  Eq.  (29)  into  Eq.  (30)  and  solving  for  dhnc 
give 

draW]  mv  dmv 

d^nc  =  (^w,  “  ~  ~  ~  dhv  —  (hv  —  hWn)  — : -  (31) 

mnr.  u  rn_. 


(29) 

(30) 
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If  Pv  =  Pvs  =  f(T)  throughout  the  dehumidification  section,  Eq.  (24)  may 
be  written  as 


or 


(32) 


(33) 


By  substituting  Eqs.  (32)  and  (33)  into  Eq.  (31)  and  defining 


dm , 

dF  =  -T- 


(34) 


Eq.  (31)  becomes 


dh. 


dh. 


dF  = 


hwx  hwD 


P  PVj  mnc  ^  w  1  “  ^W|) 


Equation  (35)  may  be  simplified  by  introducing 
hw  2  “  p  =  TW2  —  T 

and 

hv  —  h  w  jj  =  hfg 


Note:  cPw  =  1 


(36) 

(37) 


Equation  (35)  in  integrated  form  is 


T 

_  f 

1 

dhnc 

pv. 

m  y 

dhv 

J 

t2 

TW2-T 

dT  + 

P  “  Pf,  ^mnc 

)  (TW2-T) 

dT 

hfg 

HI  y 

J  (*:•  'll 

dT 

T  - 
*  w  2 

T'  mn  c 

dT  \P-PVJJ 

(38) 

Equation  (38)  may  be  solved  by  using  hnc»  pvs*  hv,  and  hfg  each  as  a 
function  of  T. 


SECTION  IV 

ANALYTICAL  MODEL  APPLICATION 


The  equilibrium  and  efficiency  equations  for  the  spray  cooler  were 
programmed  in  FORTRAN  IV  language  for  computer  solution.  Com¬ 
bustion  equilibrium  calculations  were  separately  performed  to  establish 
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the, cooler  inlet  states  using  the  method  of  Ref.  3.  These  included  the 
molecular  weights,  species  mass  fractions,  volume  of  the  nonconden¬ 
sable  gas,  and  the  gas  temperatures  based  on  certain  fuel /air  ratios  of 
conventional  hydrocarbon  turbojet  fuels. 

An  exploratory  series  of  computer  runs  was  made  to  observe  the 
effect  of  varying  one  parameter  at  a  time  in  the  equations.  In  this  way, 
results  can  be  observed  in  a  qualitative  manner.  A  realistic  grid  was 
selected  covering  a  range  of  experimental  data  for  both  the  independent 
and  fixed  inlet  conditions.  The  range  of  the  assigned  values  is  indi¬ 
cated  as  Cases  I  through  V  in  Table  I  (Appendix  II) 


4.1  PROGRAM  OUTPUT  AND  TYPICAL  RESULTS 

The  computer  program  generated  printouts  of  selected  parameters 
and  graphs.  These  parameters  included  T 3,  P3,  PV3/P3,  VnCg, 
mvi/mnc,  and  Q/mnc. 

Preliminary  computer  runs  indicated  some  difficulties  with  allow¬ 
ing  the  calculations  to  proceed  to  the  temperature  T3  =  Twl.  Solution 
of  Eq.  (17)  presents  the  possibility  of  a  negative  discriminant.  Addi¬ 
tional  discussion  of  this  negative  discriminant  and  the  required  con¬ 
straint  on  the  program  is  presented  in  Appendix  IV. 

The  results  of  the  computer  calculations  are  shown  in  Figs.  7 
through  11. 

As  the  temperature  of  the  cooling  water  is  decreased,  the  tempera¬ 
ture  of  the  exhaust  gas  at  the  exit  of  the  cooler  also  decreases  for  a 
fixed  amount  of  water  injected  (see  Fig.  7).  In  addition,  the  pressure 
drop  across  the  cooler  decreases  as  the  cooling  water  temperature  de¬ 
creases.  Therefore,  within  the  constraints  of  the  model,  where  high 
temperature  exhaust  gas  may  be  generated  and  the  pressure  drop 
across  the  cooler  may  be  critical,  testing  should  be  conducted  using 
the  lowest  temperature  cooling  water  possible. 

The  second  variable  investigated  was  the  influence  of  the  cooling 
water  injection  velocity  on  the  exit  temperature  and  pressure  (Fig.  8). 
The  injected  water  velocity  has  negligible  effect  on  the  final  gas  tem¬ 
perature  for  a  fixed  amount  of  cooling  water  injected  since  the  kinetic 
energy  of  this  item  is  almost  negligible  itself  especially  at  the  lower 
values  of  mw^/mnc.  Increasing  the  injection  velocity  will  decrease 

the  pressure  drop.  This  would  be  expected  since  the  cooling  water 
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possesses  more  momentum  and  requires  less  contribution  from  the  gas 
stream  for  its  acceleration  to  equilibrium  conditions. 

The  effect  of  decreasing  the  cooler  inlet  pressure  while  the  remain¬ 
ing  inlet  conditions  remain  fixed  is  shown  in  Fig.  9.  The  pronounced 
effect  on  the  exit  temperature  at  initial  saturation  conditions  is  caused 
by  the  low  saturation  temperature  at  the  lower  values  of  partial  pres¬ 
sure.  For  the  initial  conditions  used,  the  cooling  process  is  essen¬ 
tially  constant  pressure,  and  lowering  the  inlet  pressure  does  not 
appreciably  change  this. 

The  effect  of  increasing  the  inlet  gas  velocity  is  shown  in  Fig.  10. 
The  general  shape  of  the  exit  temperature  and  pressure  curve  is  similar 
to  that  of  Fig.  8  where  the  velocity  of  the  injected  cooling  water  is 
varied.  In  both  cases,  the  exit  temperature  and  pressure  decrease  as 
the  velocity  difference  between  exhaust  gas  and  inlet  cooling  water  in¬ 
crease.  It  can  also  be  seen  in  Fig.  10  that  as  the  velocity  increases 
there  is  a  maximum  amount  of  cooling  water  (mw^  2^mnc^  that  can  be 
injected  and  still  have  equilibrium  conditions  at  the  exit.  This  also 
places  a  minimum  value  on  the  exit  temperature  of  the  exhaust  gas. 

The  physical  significance  of  this  limit  is  discussed  in  Appendix  III, 

The  final  variable  used  in  the  study  was  inlet  exhaust  gas  tempera¬ 
ture.  This  was  varied  from  1048  to  4057°R*  and  the  exit  temperature 
and  pressure  were  determined  at  various  values  of  cooling  water  flow 
rate  (mwl  2/mnc^*  For  a  predetermined  value  of  mwl  2^mnc*  t^le 

lower  the  inlet  temperature,  the  lower  the  exit  temperature  (Fig.  11). 
The  exit  pressure  decreased  as  the  inlet  temperature  decreased,  also 
for  a  predetermined  value  of  mw^  2^mnc’  The  decrease  in  pressure 
at  a  lower  inlet  temperature  was  because  a  smaller  quantity  of  water 
was  vaporized  and  therefore  a  negligible  increase  in  stagnation  pres¬ 
sure  resulted  (A  discussion  of  the  effects  of  mass  additions  on  stagna¬ 
tion  pressure  may  be  found  in  Chapter  8  of  Ref.  4).  Also  there  was 
more  liquid  water  remaining  in  the  stream  to  accelerate  to  kinetic 
equilibrium  at  the  lower  values  of  temperature. 


4.2  KINETIC  EQUILIBRIUM 

An  assumption  made  when  deriving  the  equations  for  the  spray 
cooler  performance  program  was  that  the  flow  is  in  kinetic  and  thermo¬ 
dynamic  equilibrium  at  the  cooler  exit.  To  determine  the  effect  of 
assuming  kinetic  equilibrium,  atypical  set  of  cooler  conditions  was 


12 


AEDC-TR-71-60 


chosen,  and  the  usual  parameters  T3,  P3,  and  “wn^^nc  were  ob¬ 
tained,  as  well  as  the  six  individual  terms  of  the  momentum  equation 
in  the  following  form: 


P  lA  x  /mvl  \  ^  w  j 

— -  f  I--  T  1  -  +  - - 

mnc  \mnc  f  6c  mnc  6c 


P3  A 


where 


P  lA  l 


p  lA  l 


3  a3 


Pi 


(—A 

\  mnc  / 


^1103  mW3  ^ncj 

-  +  -  - 


mnc 


P1  Rnc  T1 


and  similarly 


'nc  Pnc  ^  ^"nc^Ai  Pnci^nci  Pnc^^ncj^ 


P3A3  P3  Rnc  T3 


P  V 

r  nc3  *  nC3 


The  assumption  of  kinetic  equilibrium  requires  that  the  liquid  water 
must  be  accelerated  to  the  gas  velocity  leaving  the  cooler.  Therefore, 
the  relative  weight  of  the  water  term  (mW3/mnc)( VnCg/g)  was  observed. 
The  inlet  parameters  were: 


and 


Pj  =  6  psia 

Ti  =  2718°R 
TWi=  520°  R 
Vw  j  =  100  ft/ sec 
Vgj  =  480  ft/sec 


mnc  =  1  lbm/sec 


Table  II  contains  values  of  the  individual  terms  which  comprise  the 
momentum  equation  and  their  respective  values  as  mw-^  2^“1nc  varied. 

The  right  hand  column  gives  the  contribution  of  the  term 
(mW3/mnc)(Vnc3/g)  to  the  total  momentum  as  the  cooling  water  ratio 

(mwl  2/]™1nc)  varied.  Since  assuming  kinetic  equilibrium  (i.  e., 

VW3  -  V11C3)  is  the  limiting  case,  the  percentages  shown  are  the  maxi¬ 
mum  contributions  of  the  liquid  water  at  the  exit  when  assumed  to  be 
in  kinetic  equilibrium.  If  the  value  of  the  water  term  were  less,  the 
effect  would  be  smaller. 


Experience  to  date  has  shown  that  a  value  of  approximately  five  or 
less  for  the  ratio  mwi+2^nc  is  sufficient  for  the  majority  of  turbojet 
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testing  conditions.  Thus,  for  the  case  presented  in  Table  II,  the  term 
directly  associated  with  kinetic  equilibrium  is  seen  to  carry  approxi¬ 
mately  6  percent  of  the  total  momentum.  It  is  concluded  that  the 
assumption  of  kinetic  equilibrium  at  the  spray  cooler  exit  is  a  justifi¬ 
able  one  for  many  cases  and  that  the  prediction  of  the  exit  pressure 
would  be  incorrect  by  no  more  than  6  percent.  This  amount  of  error 
would  probably  be  higher  for  low  pressure  or  high  subsonic  conditions. 


4.3  THERMODYNAMIC  EQUILIBRIUM 

A  study  was  also  made  to  check  the  effect  of  assuming  thermo¬ 
dynamic  equilibrium.  Assuming  equilibrium  in  this  manner  means  that 
the  temperatures  of  all  constituents  (noncondensable  exhaust  gas,  water 
vapor,  and  liquid  water)  are  equal  as  well  as  the  partial  pressures  of 
liquid  water  and  vapor  as  they  pass  the  exit  plane  of  the  cooler.  Since 
the  final  temperature  of  the  liquid  water  (Twg)  is  a  function  of  many 
items  outside  the  scope  of  this  computer  program  (such  as  nonequilib¬ 
rium  effects),  this  parameter  was  arbitrarily  varied  to  determine  the 
effect  on  the  required. water  for  cooling  (mw^  g/^nc)-  The  tempera¬ 
ture  of  the  water  at  station  3  (Twg)  was  taken  to  be: 

1«3  =  —  (TnCjj  —  TWl)  X 

where  X  was  equal  to  0,  0.  1,  0.  5,  and  0.  7.  The  limiting  case  is  X  =  0, 
which  is  the  assumption  of  thermodynamic  equilibrium.  The  partial 
pressure  of  the  water  vapor  was  taken  to  be  a  function  of  TnCg.  The 
amount  of  water  required  to  cool  the  engine  exhaust  gas  is  extremely 
sensitive  to  the  final  water  temperature  (T^)  as  shown  in  Fig.  12.  The 
input  conditions  were  the  same  as  for  the  previous  kinetic  equilibrium 
case.  To  cool  the  exhaust  gas  to  590°R  would  require  the  mass  ratio  of 
water  to  noncondensable  gas  to  vary  from  approximately  4  to  12  for  X 
varying  from  0  to  0.  7,  respectively.  Hence,  the  cooler  performance 
is  very  sensitive  to  the  degree  of  thermodynamic  nonequilibrium  at  the 
exit.  An  accurate  independent  measurement  of  the  water  and  exhaust 
gas  temperatures  is  essential  for  determining  the  true  state  of  the 
exhaust  gas  at  the  exit. 


4.4  EFFICIENCY  OF  AN  EXPERIMENTAL  COOLER 

To  show  the  effects  on  the  theoretical  efficiency, 

^theo  =  .  w  °P*imurn - ^  0f  varying  the  parameters  TWi  and  VW1,  the 

mw  equilibrium 

limiting  values  were  selected  from  Case  I  and  the  limiting  and  mean 
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values  from  Case  II,  respectively.  In  Fig.  13a,  Tw^  has  the  values  of 
500  and  545°R  which  experience  has  shbwn  are  the  lower  and  upper 
limits  of  the  cooling  water  inlet  temperature  in  the  AEDC  coolers  of 
interest.  As  would  be  expected,  the  lower  water  temperature  yields  a 
higher  efficiency  for  a  given  volume  flow  rate.  Figure  13b  shows  how 
Vw^  can  affect  the  efficiency  ;  vwl  took  the  values  50,  90,  and  130  ft/sec. 
For  a  given  volume  flow  rate,  the  efficiency  is  greatest  at  the  highest 
water  velocity.  This  is  to  be  expected.  High  water  velocity  means 
less  momentum  lost  by  the  exhaust  gas  and,  therefore,  a  situation 
closer  to  the  "ideal"  as  defined  by  constant  pressure  throughout  the 
cooler. 

The  theoretical  efficiency  (fyheo)  establishes  a  ratio  of  optimum  to 
equilibrium  conditions.  The  actual  efficiency  (hactUal^  establishes  a 
ratio  of  ideal  to  actual  conditions.  The  efficiency  ratio  (Hactual/htheo) 
establishes  a  ratio  of  actual  to  theoretical  efficiencies  resulting  in  a 
measure  of  attainment  of  the  actual  case  to  the  equilibrium  case.  This 
provides  comparison  of  actual  test  conditions  with  the  ideal  equilibrium 
conditions  and,  as  such,  is  a  measure  of  effectiveness,  of  actual  opera¬ 
tion. 


To  develop  this  efficiency  ratio  ( ^actual /’'Jtheo^  three  experimental 
tests  were  selected  on  the  basis  of  having  approximately  constant  cooler 
inlet  conditions  (see  Table  III)  and  assumed  equilibrium  outlet  conditions. 
With  the  latter  assumption  and  measuring  the  necessary  inlet  conditions 
(noncondensable  and  cooling  water  flow  rates)  and  exit  conditions  (pres¬ 
sure  and  temperature),  the  volumetric  flow  rate  of  the  exhaust  gas  at 
the  exit  is  calculated.  The  temperature  is  measured  with  a  bare  wire 
thermocouple  extending  into  the  gas  stream  where  it  will  record  the 
temperature  of  the  liquid  water  impinging  on  it.  Since  equilibrium  has 
been  assumed,  this  will  also  be  the  gas  temperature.  The  amount  of 
cooling  water  necessary  to  give  this  volumetric  flow  rate  is  calculated 
using  the  equilibrium  and  "optimum"  cooler  computer  programs.  The 
theoretical  efficiency  (ntheo  =  ^optimum^  ^equilibrium*  as  a  func_ 
tion  of  the  cooling  water  ratio  (mw^^/^nc)  is  shown  for  each- experi¬ 
mental  measurement  separately  (Fig.  13c).  This  shows  the  maximum 
achievable  efficiency  based  on  the  equilibrium  assumptions  and  the 
actual  efficiency  based  on  measured  values.  The  ratio  of  actual  to  theo¬ 
retical  efficiencies  is  shown  in  Fig.  13d  to  indicate  how  close  to  equilib¬ 
rium  efficiency  the  existing  coolers  are  operating.  It  can  be  seen  that 
the  cooler  operates  close  to  maximum  predicted  efficiency  at  the  lower 
values  of  mw^+2/mnc  but  decreases  as  the  value  of  rhwj^/^nc  hi- 
creases. 
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SECTION  V 

CONCLUDING  REMARKS 


Testing  of  turbojet  engines  and  rocket  motors  at  the  ETF  requires 
the  cooling  of  hot  exhaust  gases,  and  direct  contact  water  sprays  are 
used  for  cooling.  To  facilitate  the  understanding  of  the  cooling  process, 
an  analytical  model  of  an  exhaust  gas  spray  cooler  was  developed. 
Results  from  the  model  were  used  to  compare  with  measurements  made 
during  actual  engine  tests .  The  analytical  model  as  developed  includes 
the  equilibrium  equations  and  "ideal"  efficiency  equations  and  is  based 
on  the  assumptions  of:  (1)  thermodynamic  and  kinetic  equilibrium  at 
the  exit,  (2)  perfect  gases  and  vapor,  (3)  all  cooling  water  remaining 
in  the  stream  flow,  (4)  one -dimensional  constant  area  flow,  (5)  no  ex¬ 
ternal  drag  forces,  and  (6)  discrete  saturation  and  dehumidification 
processes.  The  method  of  solution  of  these  equations  has  been  outlined 
in  the  form  in  which  they  were  programmed  for  the  computer. 

The  computer  program  gives  an  opportunity  to  evaluate  the  influ¬ 
ence  of  various  parameters  on  the  output.  The  analytical  examination 
of  the  kinetic  and  thermal  equilibrium  assumptions  show  that  the  calcu¬ 
lation  of  cooler  exit  conditions  is  relatively  insensitive  to  kinetic 
equilibrium,  but  sensitive  to  the  degree  of  thermal  nonequilibrium. 
Efficiency  calculations  confirm  that  two  contributing  items  to  higher 
efficiency  are  low  water  temperature  and  high  water  velocity.  Experi¬ 
mental  measurements  indicate  that  the  efficiency  is  closer  to  theo¬ 
retical  at  low  values  of  mwl+  2/mnc  than  at  the  higher  values. 

o 

Verification  of  the  validity  of  the  analytical  model  assumptions  re¬ 
garding  exit  equilibrium  conditions  will  require  development  of  new 
instrumentation  for  two-phase  flow  conditions.  Further  refinement  of 
the  analytical  model  will  be  required  to  incorporate  consideration  of 
droplet  size  and  distribution,  and  nonequilibrium  thermodynamic  effects. 
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Fig.  1  Exhaust  Gas  Cooler  in  Propulsion  Development  Test  Cell  (T-1) 
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a.  Section  View  Showing  Internal  Configuration 
Fig.  2  Arrangement  of  Spray  Banks  in  T-1  Cooler 


to 

CD 


b.  Internal  View  Showing  Bank  Configuration 
Fig.  2  Continued 
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c.  Multinozzle  Spray  Heads,  View  Looking  Upstream  (Banks  No.  4  through  11) 

Fig.  2  Concluded 
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Fig.  3  System  Schematic  Showing  Instrumentation  Location 
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Water  Water 


Fig.  4  Schematic  Diagram  of  Cooler 


Fig.  5  Analytical  Model  of  Evaporative  Cooler  (Schematic) 


Fig.  6  Analytical  Model  Element  of  Dehumidification  Cooler 
(Schematic) 
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Fig.  7  Effect  of  Cooling  Water  Temperature  on  Exhaust  Gas 
Temperature  and  Pressure  (Case  I) 
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Fig.  8  Influence  of  Cooling  Water  Injection  Velocity  on 
Exhaust  Gas  Temperature  and  Pressure  (Case  II) 


34 


Pressure,  P-,  \b1 Inf  Exit  Gas 


AEDC-TR-71-60 


Fig.  9  Change  in  Exhaust  Gas  Temperature  and  Pressure  Caused  by  Varying 
the  Cooler  Inlet  Pressure  (Case  III) 
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Fig.  10  Influence  of  Inlet  Gas  Velocity  on  the  Exhaust  Gas  Temperature 
and  Pressure  (Case  IV} 


Fig.  11  Effect  of  Inlet  Gas  Temperature  Change  on  the  Exhaust  Gas 
Temperature  and  Pressure  (Case  V) 
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Fig.  12  Influence  of  Thermodynamic  Nonequilibrium  on  the  Exhaust 
Exit  Temperature  and  Pressure 


Theoretical  Spray  Cooler  Efficiency, 


a.  Effects  of  Changing  the  Cooling  Water  Temperature 
Fig.  13  Influence  of  Cooling  Water  on  Efficiency 
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Theoretical  Spray  Cooler  Efficiency, 


b.  Effects  of  Changing  the  Cooling  Water  Injection  Velocity 
Fig.  13  Continued 
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d.  Effect  of  Varying  the  Mass  Flow  Rate  Ratio  on  the  Degree 
of  Theoretical  Efficiency  Achieved 

Fig.  13  Concluded 


Parameters 


T 

l 

°R 

P 

l 

psia 

V 

nc 

i 

ft/sec 

T 

w 

°i'a 

V 

w 

1,2 

ft/sec 

2718 

5 

480 

500  to 
580 

100 

2718 

5 

480 

520 

50  to 
130 

2718 

2  to 
14 

480 

520 

100 

2718 

5 

320  to 
2130  . 

520 

100 

1048  to 
4057 


480 


520 


100 
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(COMPUTER  CASES) 


Remarks 


was  varied  in  5-deg  increments 

1,s 


V  was  varied  in  10-ft/sec  increments 
w 

1,2 


p  was  varied  in  2-psi  increments 


V  was  varied  in  nonuniform  increments 
nc 


T  was  varied  in  nonuniform  increments 


KINETIC  EQUILIBRIUM  INFLUENCE 


m 

ra 

IB 

m  V 

w  w 

1  +  2  1 

P  A 

3  3 

m 

nc 

ra 

nc 

B 

m  g 

nc  & 

™nc 

1.15 

316.4 

15.5 

3.5 

323.0 

5.13 

316.4 

15.5 

15.9 

319.8 

10.46 

316.4 

15.5 

32.4 

322.8 

23.41 

316.4 

15.5 

72.6 

332.1 

TABLE  III 

SPRAY  COOLER  CONDITIONS  USED  FOR  CALCULATING  EFFICIENCIES 


Run  No. 


6-6 

6-5 

6-23 


Pi 

psia 


Ti 

°R 

V, 

ft /sec- 

T 

Wi 

°R 

V 

f  t/sec 

m  /m 

Wi+2  nc 

Actual 

1096 

368 

512 

70 

1.42 

1096 

360 

512 

80 

2.60 

1096 

374 

512 

100 

3.16 
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APPENDIX  III 

PROCEDURE  FOR  SOLUTION  OF  EQUATIONS 


I.  Solution  of  Equilibrium  Equations 

The  procedure  for  determining  exit  conditions  by  solving  Eqs.  (7) 
or  (11)  and  (17)  involves  two  separate  steps  as  follows: 

For  that  part  of  the  process  leading  to  saturation: 

1.  Let  P2  =  Pi. 

2.  Assume  a  T2. 

3.  Calculate  a  new  temperature  (T2new)  using  Eq.  (7), 

Eqs.  (9)  and  (10)  in  which  Dalton's  law  is  used  to 
eliminate  Pnc>  and  the  partial  pres  sure -temperature 
relationship. 

4.  Compare  T2new  with  the  assumed  temperature  in 
step  2.  If  their  difference  is  within  tolerance,  pro¬ 
ceed  to  step  5;  otherwise  let  T2  =  T2new,  and  go  back 
to  step  3. 

5.  Calculate  a  new  pressure  (p2new)  using  Eq.  (17). 

6.  Compare  p2new  with  p2  in  step  1.  If  their  difference 
is  within  tolerance,  P2  is  the  saturation  pressure,  and 
T2  is  the  saturation  temperature;  otherwise  let 

p2  =  P2neW'  T2  =  T2new.  and  go  back  to  step  3. 

For  that  part  of  the  process  beyond  saturation  (dehumidification): 

1.  Let  T3  =  Tsaturatjon  -  AT. where  AT  is  some  small 
temperature  increment  on  the  order  of  one  degree. 

2.  Let  P3  =  Ps. 

3.  Calculate  a  new  pressure  (p3new)  using  the  partial 

pressure-temperature  relationship,  Eqs.  (9),  (11), 
and  (17). 

4.  Compare  p3new  with  the  pressure  (P3)  in  step  2.  If 
their  difference  is  within  tolerance,  the  equations 
are  satisfied  at  the  assumed  temperature  (T2).  Go  to 

p3  +  p3new 

step  5.  Otherwise  let  P3  = - 2 - ,  and  go  to 

step  3. 
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5.  Let  T g  =  T 3  _  AT.  Go  to  step  2. 

The  dehumidification  procedure  is  allowed  to  continue  until  Tg  =  TWg* 


II.  Solution  of  Efficiency  Equations 


The  procedure  for  solution  of  Eq.  (27)  for  the  saturation  case  is  as 
follows : 


1. 

2. 

3. 

4. 


5. 


Determine  PVg,  a  known  function  of  T,  at  an  assumed  T3. 

Calculate  Pvq  from  Eq.  (28). 
jnew 


Determine.  T3new,  a  known  function  of  Pv>  at  Pv3new* 

Compare  Pv3new  with  PVg.  If  their  difference  is  within 

tolerance,  let  PVo  =  Pv_,  and  proceed  to  step  5;  other- 
3new  v3 


P  +  P 
v3  v3 


new 


wise  let  PVo  = - - -  and  To  =  To  ,  and  go 

v3  2  >5  ■mew 

back  to  step  2. 

Calculate  the  water  required  for  stauration  using 
Eq'.  (27). 


For  the  dehumidification  case,  the  solution  of  Eq.  (38)  is  required. 
The  procedure  for  solving  this  equation  is  based  on  a  Gaussian  integra¬ 
tion  technique.  The  term  (F)  is  the  ratio  of  the  mass  flow  of  the  water 
injected  into  the  dehumidification  cooler  element  to  the  mass  flow  of 
noncondensable  gases,  F  =  mw^/^nc-  To  determine  the  total  mass 
flow  ratio  of  the  entire  cooler,  the  water  injected  into  the  saturation 
section  (mW2evap/ *nnc)  can  be  added  to  F.  This  will  give  the  ideal 

mass  flow  ratio  (niw^g^/mnc).  From  this  can  be  calculated  the  value 
of  total  water  required  for  the  entire  cooler  for  the  idealized  case 
(^widead.)- 
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APPENDIX  IV 

SOLUTIONS  OF  THE  MOMENTUM  EQUATION 


The  original  development  of  the  momentum  equation  in  Section  III 
led  to  the  solution  of  the  exit  pressure  (P3)  from  an  equation  of  the 
quadratic  type: 


where  a,  b,  and  c  have  previously  been  defined.  The  possibility  of  the 
quantity  b2  -  4ac  becoming  negative  and  giving  an  exit  pressure  in  terms 
of  a  complex  number  is  immediately  recognized.  The  question  becomes 
( 1)  will  this  happen  during  numerical  solution,  and  (2)  if  it  does,  is 
there  a  physical  explanation? 

For  certain  cooler  inlet  conditions,  b2  -  4ac  does  become  negative 
during  the  numerical  solution.  Correspondingly,  there  is  a  minimum 
value  of  T2  and  a  maximum  value  of  mw^  2^nc»  say  (mW]_  2^”lnc^max* 

for  which  a  solution  can  be  obtained  to  the  momentum  equation.  A 
physical  interpretation  of  this  behavior  is  that  it  is  impossible  for 
equilibrium  to  exist  at  the  exit  of  the  cooler  when  the  value  of  the  water 
added  to  the  cooler  exceeds  (mwlj  2^nc^max-  the  cooler  is  oper¬ 
ating  with  equilibrium  at  the  exit  and  the  rate  of  cooling  water  injection 
equals  to  (mwlj  2/r”nc)max  and  the  rate  of  cooling  water  injection  is  in¬ 
creased,  the  gas  phase  will  choke  at  the  cooler  exit,  and  thermodynamic 
and  kinetic  equilibrium  cannot  exist  at  the  exit.  Any  further  increase  in 
water  injection  rate  must  cause  an  adjustment  in  the  initial  conditions. 
The  analysis  indicates  a  sudden  choking  condition.  That  is,  the  gas 
phase  Mach  number  at  the  cooler  exit  would  suddenly  jump  from  some 
subsonic  valve  to  unity  with  the  addition  of  a  differential  amount  of  water. 
However,  in  an  actual  cooler  which  has  finite  length,  sudden  choking 
would  not  be  expected  to  occur  because  of  nonequilibrium  effects.  Also, 
the  value  of  mwi>  2^nc  a*  which  choking  would  occur  in  an  actual  cooler 
would  be  somewhat  higher  than  (mw^  2^r^1nc^max- 
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